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Abstract. In this paper, a new formalism for the filtered density function (FDF) approach 
is developed for the treatment of turbulent polydispersed two-phase flows in LES simulations. 
Contrary to the FDF used for turbulent reactive single-phase flows, the present formalislm 
is based on Lagrangian quantities and, in particular, on the Lagrangian filtered mass density 
function (LFMDF) as the central concept. This framework allows modeling and simulation of 
particle flows for LES to be set in a rigorous context and various links with other approaches 
to be made. In particular, the relation between LES for particle simulations of single-phase 
flows and Smoothed Particle Hydrodynamics (SPH) is put forward. Then, the discussion and 
derivation of possible subgrid stochastic models used for Lagrangian models in two-phase flows 
can set in a clear probabilistic equivalence with the corresponding LFMDF. 
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1. Introduction 

Dispersion and entrainment of inertial heavy particles in turbulent flows are crucial in a number 
of industrial applications and environmental phenomena. Examples are mixing, combustion, 
depulveration, spray dynamics, pollutant dispersion or cloud dynamics. Direct numerical 
simulation (DNS) together with Lagrangian particle tracking (LPT) has been successfully used 
to investigate and quantify the behavior of particles in turbulent flow. DNS studies have shown, 
for instance, that particle segregation and subsequent deposition phenomen a are related with the 
inter action between particles and the near- wall coherent vortical structures ( Soldati &: Marchiolil . 
2009). Furthermore, large-scale clustering and preferential concentration, due to the tendency 
of inertial particles to distribute pr eferentially at th e periphery of strong v o rtical regions and 
to se gregate into straining regions (|Bec et al\ . 120071 : IWang fc Maxevl . 1 19931 : iRouson fc Eatonl . 
200 ll ). have been revealed. 

On the other hand, DNS remains limited to low Reynolds numbers and simple geometries. 
Conversely, thanks to the huge growth of available computational resources, large-eddy 
simulation (LES) can address increasingly complex problems. In LES the basic idea is to 
directly simulate only the turbulence scales larger than a given dimensi on, while the e ffects of the 
unresolved subgrid scales (SGS) on the large-scale motion are modeled( Sagaut . 20061 ). Thus, the 
exact dynamics of the filtered fluid velocity field can at best be obtained from LE S (ideal LES) , 
this when SGS mode ling and numerical errors are negligible. Recent studies ( Kuerten . 20061 : 
Marchioli et al . 20081 ) have demonstrated that neglecting the effects of the fluid SGS velocity 



fluctuations on particle motion leads to significant underestimation of preferential concentration 
and, consequently, to lower deposition fluxes and reduced near-wall accumulation. The issue 
is therefore how to model the SGS effects in the particle equations. Actually, this i ssue is 
very similar to the one appearing in the modeling of single-phase turbulent reactive flows (jPopd . 
1994) and, in this case, stochastic models that can treat the reactive source terms without 
approximation have shown their great potential. For the same reason, a stochastic approach 
to polydispersed turbulent two-phase flows is interesting. First ideas of stochastic models have 
been proposed for LPT coupled with LES of the carrier phase, but they ha ve mainly been used 
for homogeneous isotropic turbulence ( Pozorski Aptel . 20091 ; Gobert . 20 id ) and they have been 
proposed with no clear theoretical background. 

In this paper, we develop a rigorous probabilistic formalism within which one can 
develop s tochastic models fo r turbulent two-phase flows. As put forward by Pope and 
coworkers) Sheikhi et al\ . 20071 ). this modelling requires the definition of a proper Filtered density 
function (FDF). However, at variance with Pope's case, we develop a Lagrangian formalism 
which includes naturally the Eulerian one. In this way, we also show how LES is related to 
another Lagrangian method, smoothed particle hydrodynamics. 



2. Governing equations 

We consider in this work heavy particles carried by a fluid. The equations of motion of the fluid 
are Navier-Stokes equations 



dUj v dUj 1 dP ^ d 2 Uj 
dt 3 dxj p dxi dx 



+ U ^ = -^ + ^ ( lb ) 



iF +t, '&r r W+ s <<*> (lc) 

where the set of scalars 3? = (mi, . . . , mjv s , T) are the mass fractions mi(l = 1, . . . , N s ) of the 
N s species that compose the reactive mixture and the temperature T (if we consider turbulent 
reactive single-phase flows). For heavy particles where p p 3> pj transported by the fluid flow, 
the drag and gravity forces are the dominant forces and the particle equation of motion can be 
simplified to 

^ = i(U,-U p)+g . (2) 

The drag force has been written in this form to bring out the particle relaxation time scale 

_ _ p p 4 d p 
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(3) 



The drag coefficient is an empirical coeffi cient that can b e estimated through experiments. 
Various expressions have been put forward( Clift et al . 19781 ). such as 
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[ 1 + 0.15i?e°- 687 ] if Re v < 1000, 



C D = I Re p " P j , (4) 



0.44 if Re v > 1000. 



3. Eulerian and Lagrangian descriptions in PDF methods 



In the case of compressible or reactive flows, the proper probability density function (pdf) 
to be considered is the mass density function Fl rather than the pdf itself (which is instead 
the right choice in the incomp r essibl e case since all fluid particles have a constant mass for 



a given volume) (jSheikhi et all [2007) . Since particles behave like a compresible fluid and 
may display polydisper s ion, t he mass density function is also a relevant choice for two-phase 
flows ( Minier Peirano . 200 ll ). In this section, we r e call s ome basic definitions (a thorough 



description can be found elsewhere) Minier Peirano . l200ll )). For a complete description, we 



introduce the mass density function F^t; y, V, i Sf) where F^it; y, V, \l>) dy dV d^Sf represents the 
probable mass of fluid particles contained in an element of volume dy dV d& in phase space. 
The mass density function is consequently normalized by the total mass of fluid 

M = j F L (t- y, V, *) dy dV d*. (5) 

These relations are valid both for fluid and solid particles. The correspondence between the 
Eulerian and the Lagrangian descriptions is given by, 

F £ (t,x;V,*)=F L (t;y = x,V,*) = J F L (t; y, V, *) <5(x - y) dy, (6) 

where Fe is the Eulerian mass density function. This relation simply expresses the fact that, 
when in phase space y = x, there is equivalence between the two points of view since we 
consider the same event (in the incompressible case this relation is verified by pE(t,x; V) and 
PL(t; y, V)). Mass continuity constraint imposes that the integral of Fe over phase space (V, 
is the expected density at (t, x) (the probable mass of particles in a given state per unit volume). 
The expected density, denoted (p)(£,x), is thus defined by 

a(i,x)(p)(i,x) = J F E (t,x;V,y)dVdV, (7) 
3.1. PDF theory: No-model equation 

We recall the exact equations for the trajectories of fluid and discrete particles in a Lagrangian 
form. The new set of equations take the following form, 

'/•'•/, 'V.;'/'- 

dU+. = A+.dt + A+^ f (t,Z,(Z))dt, 
d*+ =I>A*+ dt + S f (*+)dt, 
dr + ■ = U + - dt 



du; t i = A i,i dt > 



d^;, k = r p A^ k dt + s p (^)dt, 

where the indexes I and k refer to the dimensions of <&/ and <& p , respectively. The + subscript 
is used to indicate the exact trajectories in contrast to the modelled ones. Both (exact) 
accelerations are given by (Aj i and A^ i are the accelerations of the fluid particles and the 
discrete particles, respectively) 



U ~ 7f~d^ + U f,i ' (9) 
^ = ^+-11+) + 9i . 



A new term is added in the momentum equation of the fluid to account for the influence of the 
particles on the fluid. 

Using the techniques developed in stochastic calculus, it is possible to show that the transport 
equation verified by the Eulerian mass density functions Fe for both phases is of the form (here 
we show only that for solid particles for the sake of simplcity): 

g F P apP a 



8t **dy Vii dv t 



° {[A Sii + {Ap^i] y p , V p , *„)] F% ) (11) 
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4. The LES approach for single-phase turbulent flows 

In single-phase flows, the starting point is provided by the Navier-Stokes equations supplemented 
with conservation equations for a set of scalars. These are field equations (the different fields 
are density p(t,x), pressure P(t,x), velocity U(i, x) and scalars <fr(t,x), where x represents the 
coordinates in physical space). The LES approach consists in calculating the larger three- 
dimensional unsteady turbulent motions while the effects of the smaller scale motions are 
modelled. For that purpose, LES involves the use of a spatial filter 

/oo 
/(y,t)G(y,x)dy (12) 
-oo 

where G is the filter function, f(x,t)L represents the filtered value of the transport variable 
f(x,t), and /' = / — denotes the fluctuations of f from the filtered value. We consider 
spatially and temporally invariant and localized filter functions, thus G(y, x) = G(x — y) with 
the properties, G(x) = G(—x), and J G(x)dx = 1. 

Starting from Navier-Stokes equations, application of the filtering operator yields: 



dxj 



(13) 



d(Ui)L , r v d(Uj)L _ ^d(P) L d( Tlj ) L 

+ {u > h ^x--—p^x- + uAmL -^x-> (14) 
+ Wi-teT " { h — + {Sl) > (15) 

where {t^l = {UiUj) L - {Ui) L {Uj) L and (M$) L = (*jL7j-) l - {^i) L {Uj) L denote the subgrid 
stress and the subgrid mass flux respectively. 

4-1- The Filtered Density Function (FDF) approach 

In this section, we recall the main points of the FDF theory introduc ed for reactive flow s, 



following the approach proposed by Pope and developed with coworkers (|Sheikhi et all 120071 ). 
For general compressible single-phase flows, the FDF theory introduces a "filtered mass density 
function" (FMDF), denoted F, 



F(t, x; V, *) = / p(U(y, t), #(y, t), V, *)G(y - x)d y , 



(16) 



p(U(y, t), *(y, t), V, *) = 5(V - U(y, t)) ® S(* - 3>(y, t)) , (17) 



where 5 denotes the delta function and denotes the composition domain of the scalar array. 
The term p(U(y,t),$(y,t) is the fine-grained mass density and Eq. (|16|) implies that the 
FMDF is the spatially filtered value of the fine-grained mass density Thus, F gives the density 
in the composition space of the fluid around x weighted by the filter G. With the condition of a 
positive filter kernel, F has all the properties of a PDF. It is also useful to define the conditional 
filtered value of a variable Q(x,t): 

(Q(x, t) i r ,. )h . JQ(y - t)p( ^ ( ;™ )G(y - x)rfy , (18) 

which implies the important integral property: 

(Q(x, t)) L = J (Q(x, t) | *) L F(t, x; V, V)dVdV. 
By standard techniques, it is possible to show that the governing equation for the fdf is 

f + W = -^ {A '- 1 v - *> F] - A 1 (rA *' 1 v - *> F1 - A ls ' w n (19) 

The important aspect is that only an Eulerian FMDF i s considered and th at, apart from the 
numerical Monte-Carlo particles used later on the theory ( Sheikhi et al . 20071 ). there is no notion 
of a Lagrangian filtered density. 



5. A new FDF formalism for polydispersed two-phase flows 

In turbulent reactive single-flows, the FDF formalism developed so far by Pope and co-workers 
is based only on the Eulerian filtered density function, which is consistent with the fact that the 
basic exact equations to handle are field equations. However, in polydispersed two-phase flows, 
the physical situation is completely different since the basic exact equations (for the particle 
phase) are Lagrangian. Therefore, in order to develop a proper LES treatment for particle 
flows, it is not possible to copy simply what has been achieved in single-phase flows. A new 
point of view must be introduced and a new formalism must be entirely developed. 

The theory briefly recalled in Section [3] has shown that the key notion is the Lagrangian 
Mass Density Function. For the development of a FDF approach to two-phase flows, we choose 
therefore to retain a Lagrangian point of view as the starting point and we introduce as the 
central concept the Lagrangian filtered mass density function (LFMDF) which is defined as 
(considering, for sake of simplicity, only a particle phase) 

Fl(t;y,V,*) = J F p L (t;y',V,*)G(y-y')dy' , (20) 

In terms of fine-grained densities, this can be written for the description of the particle phase 
(noted F p ) with iV the number of particles present in the domain at the time t as: 

r N 

F p L {t-y,V,*) = / ^m PI^ G(y-y')<5( ? /'- a;^ (^))®5(F-^7 ^ (^))®5(*-* ^ (^)) ^ ^y , 

i=l 

N 

= ^m p>i G(y-x i (t))®(5v'®53, . (21) 

i=l 



In order to avoid a confusion between the filtering operator and the Eulerian and Lagrangian 
indices (E and L, respectively) which will have to be introduced, the filtering operator of a 
variable A will be noted by A (instead of (A)l in previous sections). 

For further developments, it is useful to consider the conditional filtered value of a variable 
H(t) which is defined as 

*)Fi = j H(t)Fl(t; y', V, *)G(y - y')dy' , 

N 

= ^2 Hijn pA G(y - x ; ) ® 5y <8> • (22) 



i=l 



As in Section f4.il this operator has all the properties of a conditional average and a key property 
is that for any function of the variables in the state-vector, H = H(x, U, <&) 

(Hfcv,*)=H(y,V,*). (23) 

In this formalism, following what is done in PDF methods, we choose to define the Eulerian 
filtered mass density function (EFMDF) directly from the LFMDF as 

f*(t,x;V,*) = *£(t;y = x,V / ,*). (24) 
From these definitions, we have thus that 

**(t,x;V,*) = J F£(i;y',V,*)G(x-y')dy' . (25) 
or in terms of the fine-grained densities 

TV 

F|(i,x;V,*)= J2m P)i G(x-x 1 (t))®5v®5*. (26) 

i=l 

From the EFMDF, we can define the filtered for any variable Q(t,x) 

a p (t,x)(p) p Q(t,x) = J J (Q\y~y,*) F*(t,x; V, ^)dVd^ (27) 
and for a function of the variables in the state-vector 

Op(t,x)(p)pQ(t,x) = J j Q(t,x,V,V) ^(t,x;V,^)dVcty (28) 

where a p (t,x){p) p is the filtered local value of the particle mass fraction at t and position x 

a k (t,y){p) p = jF%(t t x]V k ,i/> k )dV p dil> k . (29) 

Given the filtered mass density function, it is possible to derive the associa ted evolution 
equation using the fine-grained definition (f2"Tj) with standard techniques Pope ( 19851 ) (considering 



in the following derivation that m Pt i are constant, for the sake of simplicity) 

dF[ dG. My . . \ 

-gf = l^ [ m p,i-Q]7 6 v* + ™ p ,iG—5* + m p<i G—5 v j 



i=l 
N 



E( dG dx,i „d\Ji 8S V s „d<&i d5y \ 

i=l ' 

E x / dG dxi dUi d5 v „d®i d5^ \ 

= _ ±[(A^V, <t)f* ] - A[ (A,i^V, *)!*] (30) 

This equation corresponds to the system of Lagrangian equations which are written, retaining 
general notations since we will apply them for fluid and solid particles in following sections, as: 

dxi = U{ dt, (31) 
dUi = A Uti dt (32) 
d$z = dt, (33) 

Since the EFMDF is by definition identical to the LFMDF, the same equation is satisfied 
automatically by Ft, from which any transport equations for filtered fields can be easily 
extracted by integration over the corresponding sample-space, cf. Eq. (|27|) , Then, by 
decomposising into filtered and "fluctuations" (or unresolved values), the equation for the 
FMDFs can be expressed as 



-£[(v-e)3i 



^[{(A,STv.*)-A„}f»] (34) 
A[{(A t |^.*)-A;}Ff] 



where the first line corresponds to the effects of resolved scales (closed terms) and the last three 
ones to the effects of the unresolved modes (unclosed terms). 

6. Smoothed Particle Hydrodynamics (SPH) as a particle-based LES 

As a first application of the present general formalism, we consider the case of N fluid particles 
(in which case, we will use here an index / instead of p). If we disregard the effects of unresolved 
scales, the transport equation for the EFMDF is simply: 



8F[ d[UF[] _ d 



^ + U -^ = -^[A u F[] (35) 



The particle (or Monte-Carlo) solution of this equation is 



dxf t i = U f t i dt, (36) 
dU fji = Af ti dt (37) 

where Af i is the Navier-Stokes equation in a Lagrangian formulation 

A » = -7fZ + uAU »- (38) 

By applying the general relations for a particle labelled a (and using label b to express the sum 
over other particles), we have that 

UaPaA a = m b A b W ab (39) 
b 

where we have used for the filtering function, G a kernel function, W , with a typical length h 

G(y a - y b ) = W(y a - y b , h) = W ab . (40) 

For incompressible flows, the density is constant p a = p b and the local mass fraction estimation 
gives (writing m b = p b V b with V b the particle volume) 

a a Pa = ^2 m f> W ab = ^ Pt>VbW ab , and that a a = ^ V b W ab = 1 , (41) 
b b b 

which is the normalisation condition for the kernel W and states that the domain is continuously 
filled with fluid particles. Thus, we obtain the following particle equations 

<E fti = U fti dt, (42) 



dUf,i = £ 



Pf dxij b 



V b W ab dt (43) 



Therefore, by disregarding unresolved effects and considering a Lagrangian description 
of a singl e-phase flu i d flow , we obtain exactly the smooth particle hydrodynamics (SPH) 
forma lism (Monaehan. 120051 ) . Further SPH expressions, in particular for truly incompressible 
flows ( Szewc et al . 20 111 ), give the complete equations by expressing the pressure- gradient and 
viscous terms with a particle formulation: 

*- e « (44) 

dt ^ [ paPb Pa+Pb T l ah J 

Thus, as a first conclusion, it is now rigorously demonstrated that the SPH method is 
equivalent to performing a Lagrangian LES without a subgrid model. This is a rather intuitive 
notion, considering the smoothing nature of the SPH operator, but there was previously no 
theoretical framework to express it in a strict setting. In other words, present SPH formulations 
appears a particle-based LES without any sub-kernel model. And, consequently, a direct 
extension could be to consider sub-kernel models, for example with the form of particle stochastic 
models, which would provide a direct link with particle simulations in LES. 



7. LES for particle flows and subgrid stochastic models 



We now consider poly dispersed solid particles embedded in a turb ulent described by a LE S 
simulation. Following the discussion and choice made in PDF methods iMinier &: Peirano fl2Q0lh . 



it appears interesting to include the velocity of the fluid seen by the particles in the state-vector 
retained and, thus, to handle the LFMDF F^(t; y p , V p , V s ). The transport equation for the 
LFMDF (and EFMDF F%(t, y p ; V p , V s )) is then 



dt 8y p r p J L ' dV. 

which corresponds to the particle equations 

dx p = Up dt 

dVp = Us-u P dt (47) 

T p 

d\J s = A s dt. (48) 

where a model for A s has to be proposed. 

The first idea that was applied in LES simulations with solid particles was to disregard subgrid 
effects and to identify the (instantaneous) fluid velocity seen, U s (t), with the available filtered 
value at the particle location, XJ s (t) = U s (i,x p (i)). This amounts to considering only particle 
location and velocity in the corresponding LFMDF since U s (t) is now given, F p (t; y p , V p ) whose 
evolution equation is 

riytv^U^f^zaW (4 9) 



dt dy p 8V 



v \ p / 

This equation involves only first-order derivatives and is an example of a classical Liouville 
equation (note that the notation for the sample value of particle velocity, U p , has been changed 
to V p to indicate that, in this case, particle velocitues are only driven by a known drag 
force). Indeed, the Montec-Carlo solution of this transport equation corresponds simply to 
the deterministic particle equations of motion 

dxp = Up dt, (50) 

d\J p = Us I Up dt (51) 
T p 

where t p represents the particle relaxation timescale expressed only in terms of the filtered 
relative velocity U r . 

However, in order to take into account the explicit effects of unresolved scales, the next step 
is to introduce a stochastic model for U s which, in the corresponding sample-space for the 
LFMDF Ff(t; y p , V p , V s ) means that we have now a modeled transport equation, for example 
a Fokker-Planck equation 

gg | d[ V P ,Fl ] _ d ( U s , - U Pti \~ P] gj DsM } { l d 2 [ ggg ] 



dt dy Pti dV p /\ rp J LS dV s ,i 2 d 2 V s ^ 3 



which corresponds to a stochastic diffusion process for the modeled particle variables 



dx p = U p dt 

dU p = Vs ~ Up dt (53) 

Tp 

d\J s = B s dt + B s dW (54) 

with tfW the increments of the vector of independent Wiener processes. For practical examples 
of precise diffusion stochastic processes (stochastic modeling being not necessarily limited to 
this sub class!), it is possibl e to consider the extensions of the stochastic models used in PDF 
models ( Minier et al . 20041 ) . For example, a simple model would be 



dx P) i = Up^i dt (55) 

dUp . = U ^ ~ Via dt ( 56 ) 

Tp 

dU s ,i = -l^dt+(v^-UfJ)^dt- U *- dt + JcfkdWi (57) 



where crossing-trajectory effects have been kept for the timescales of the unresolved fluid velocity 
seen r£ ■ (with t£ the timescale of the fluid unresolved motions). It is thus important to note that 
the diffusion matrix, B s , is diagonal but not isotropic. In LES simulations, the fluid subgrid 
fluctuations can be regarded as nearly isotropic (as a first step) and the general expression 
proposed in PDF models for C* would be simplified to 

C* = C b i + l(bi-l) with b % = ^- (58) 

The development of the corresponding filtered particle fields, as well as a more detailed discussion 
of several possible subgrid kernels, will be presented in forthcoming manuscripts. Yet, the present 
framework already emphasizes that it is crucial to consider satisfactory diffusion coefficients (or 
diffusion matrix) for the velocity of the fluid seen. If some arguments have been put forward 
on the form of the drift vector, D s , an essential aspect remains an acceptable form of B s so as 
to reproduce, at least, a correct energy flux from the resolved to unresolved scales. Note that, 
the simple idea to use a closure expression as in single-phase flows, that is B s ^j = V Co e dij , is 
drastically false since a physically-incorrect and spurious energy flux, or transfer rate, is induced 
inconsistent with the modeled rate ?. 



8. Conclusions 

In this paper, a theoretical formalism for the LES of particle turbulent polydispersed two-phase 
flows has been presented. The formalism has been developed for the Lagrangian FDF, for which 
the general transport equation is derived. In this framework, it is possible to build stochastic 
models for particle flows in a rigorous way. Furthermore, we have applied this formalism to 
show that SPH is a particle-based LES. Finally, a first example of Langevin model constructed 
within the formalism is proposed considering isotropic sub-grid fluctuations. 
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